The Silicon Inversion Layer With A Ferromagnetic Gate: A Novel Spin Source 



J. P. McGuire, C. Ciuti, L. J. Sham 
Department of Physics, University of California, San Diego, La Jolla, CA 92093-0319. 

(Dated: February 6, 2008) 

Novel spin transport behavior is theoretically shown to result from replacing the usual metal 
(or poly-silicon) gate in a silicon field-effect transistor with a ferromagnet, separated from the 
semiconductor by an ultra-thin oxide. The spin-dependent interplay between the drift current (due 
to a source-drain bias) and the diffusion current (due to carrier leakage into the ferromagnetic gate) 
results in a rich variety of spin dependence in the current that flows through such a device. We 
examine two cases of particular interest: (1) creating a 100% spin-polarized electrical current and 
(2) creating a pure spin current without a net electrical current. A spin-valve consisting of two 
sequential ferromagnetic gates is shown to exhibit magnetoresistance dependent upon the relative 
orientations of the magnetization of the two ferromagnets. The magnetoresistance ratio grows to 
arbitrarily large values in the regime of low source-drain bias, and is limited only by the spin-flip 
time in the channel. 



The creation, manipulation, and detection of spin- 
polarized carriers are the essential ingredients of semi- 
conductor spintronicsi 1 ! 2 ! 3 Although well established us- 
ing optical methods;^ 6 - the ultimate goal is to achieve 
this level of control over the spin degree of freedom in all- 
electrical devices^ To this end, spin injection from fer- 
romagnets into semiconductors&Siifliii and spin control 
through the Rashba spin-orbit effec1ii2iiiiiii 5 rii& have at- 
tracted much interest during the last few years. In addi- 
tion, it has recently been proposed^ and experimentally 
verifiecUSiiS that a pure spin current, without an accom- 
panying electrical current, can be created in direct-gap 
semiconductors using the interference of one- and two- 
photon absorption processes. Such pure spin currents 
will allow for new probes of spin-dependent phenomena 
in semiconductors and may lead to alternate spintronics 
device designs not present in existing charge-based tech- 
nology. 

In this paper we present a method of creating fully 
spin-polarized charge currents and pure spin currents 
using a ferromagnetically-gated inversion layer in the 
regime of low source-drain biasi2Si2i Our method is pos- 
sible in any semiconductor, notably silicon, and not just 
those with well-defined optical transitions. It relies on 
a proximity effect in which semiconductor carriers ac- 
quire spin-dependent properties through coupling with 
a nearby ferromagnet^ Experiments have shown that 
optically-pumped unpolarized electrons in a semiconduc- 
tor acquire a net spin-polarization in the presence of an 
interface with a ferrromagneti2Si2i24 The magnitude of 
the spin-dependent effects are optimized by confining the 
semiconductor carriers at the interface with the ferro- 
magnet, with the ferromagnet acting as the electrically- 
biased gate responsible for the inversion. The crucial 
parameter is the thickness of the oxide barrier which 
separates the semiconductor from the ferromagnet; since 
the coupling of the two depends exponentially on their 
separation, our proposal will be most relevant for ultra- 
thin gate oxides at the forefront of current technology. 26 
The balance of growing oxides thin enough to produce 
ample coupling while preserving the functionality of the 



inversion layer is delicate, although we point out that 
we are merely taking advantage of the aggressive scal- 
ing of field-effect transistors to the nanoscale. In addi- 
tion, we will discuss the behavior of a recently-proposed 
spin- valve with two sequential ferromagnetic gates in the 
same regime of low source-drain bias, where the magne- 
toresistance ratio due to the relative orientation of the 
ferromagnet magnetizations can grow to arbitarily large 
values. 

To demonstrate the operational principle of the device, 
consider the spin-dependent current in a two-dimensional 
electron gas (2DEG). We will concentrate on the silicon 
inversion layer with an SiC>2 barrier and a ferromagnetic 
gate, and consider only a single occupied subband in the 
silicon. For simplicity, we will assume the two spin chan- 
nels are completely decoupled, but we will return to ad- 
dress any coupling between the two spin channels later in 
the paper. The device is shown in Fig.^a); the growth 
axis is in the z-direction, the in-plane current flows in the 
x-direction, and we assume the device is translationally 
invariant in the y-direction. The in-plane electrical cur- 
rent J a {x) in spin channel a = {|, J.} contains both drift 
and diffusion terms, 

d 

Ja{x) = 9a(x)E x + eD a {x)—N (x) , (1) 

where — e is the electron charge, g a (x) is the Drudc 
conductivity, E x is the electric field that drives the in- 
plane current, D a (x) is the diffusion constant, and N a (x) 
is spatially-dependent 2DEG density. To achieve spin- 
dependent transport, either the transport coefficients or 
the 2DEG density must vary considerable between the 
two spin channels. We show below how the coupling be- 
tween the 2DEG and a ferromagnetic gate achieves this 
splitting. 

In Ref. |53 the coupling was calculated from the wave- 
function matching conditions, resulting in an equation 
which specified the spin-dependent complex energy of 
the confined inversion layer state coupled to the ferro- 
magnet. In Ref. an effective tight-binding Hamilto- 
nian was derived from the effective-mass Hamiltonian of 



the silicon/Si02/ferromagnet interface. The form of this 
Hamiltonian allowed for the calculation of the 2DEG self- 
energy due to the interaction with the ferromagnet, the 
approach we will adopt in this paper. The results of 
the derivation in Ref. |2l| are summarized below in or- 
der to clarify the method of calculation used throughout 
the rest of this paper. The effective Hamiltonian of the 
2DEG coupled to the ferromagnetic gate is 



- fm a 



(2) 



A one-dimensional calculation is possible for the cou- 
pling calculation because the conservation of the in-plane 
wavevector is relaxed The operator a\ creates the 
inversion layer state with spin a, energy e sl , and wave- 
function x sl (z), which is the lowest bound state of the 
potential 

V s \z) = [-U h - eEf ■ (z - z b )] Q(z - z h ) , (3) 

where Ub is the oxide barrier height from the bottom of 
the semiconductor potential, Ef is the electric field in the 
semiconductor, and z^ is the thickness of the oxide bar- 
rier. The "0" of energy has been put at the top right side 
of the barrier to simplify the calculations. The operator 
cl creates a ferromagnet state with spin a, wavevec- 
tor k, energy ej™., and wavefunction </>j™.(.z), which are 
states of the potential 



V* m (z) = 



U im + —a ■ M 



-eE 



" (z - z h )Q{z)Q(z h - z) , (4) 



where E/f m + A/2 is the Fermi energy of the majority spin 
band in the ferromagnet, U[ m — A/2 is the Fermi energy 
of the minority spin band in the ferromagnet, and E^ is 
the electric field in the barrier. Note that for simplic- 
ity we treat the ferromagnet as exchange-split parabolic 
bandsi^ The coupling between the 2DEG and the ferro- 
magnet is approximated by 



dz [^{z)]* V^{z) X s \z) , 



(5) 



which is basically an overlap integral with an exponen- 
tially decaying term. Eq. @ is an approximation of the 
effective-mass Hamiltonian of the Si/SiO^/ferromagnet 
interface in which only the most signifcant terms in the 
coupling have been retained. 

From the simple tight-binding Hamiltonian, Eq. J5J, 
the self-energy of the 2DEG electrons due to the interac- 
tion with the ferromagnetic gate is 
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(6) 
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FIG. 1: (a) The silicon field-effect transistor with a single 
ferromagnetic gate. The source-drain bias V s d creates the in- 
plane field E x that drives the drift current. The ferromagnetic 
gate is biased (V g ) to create the inversion layer, (b) The spin- 
valve with two adjacent ferromagnetic gates; the two ferro- 
magnets have different coercive fields to enable the switching 
of the second gate's magnetization without affecting that of 
the first gate. The gap between the two ferromagnets is filled 
with a paramagnetic metal to ensure uniform 2DEG confine- 
ment throughout the device. 



To account for strong scattering in the ferromagnet, 
a large imaginary part has been added to the ferro- 
magnet states (7;[™ ). Calculating the complex energy 
E which satisfies E - e si + i0+ - E a (E) = 0, we 
can evaluate both the spin-dependent 2DEG level shift, 

A a (E) = Re[s CT (^)j, and the spin-dependent 2DEG 



lifetime, T a {E) = -ft/2Im Y,„(E) 



due to the coupling 



with the ferromagnet. The spin-splitting due to the level 
shift, |Aj — A J, is very small compared to the Fermi 
energy of the 2DEG for even the thinnest practical ox- 
ides, and will be neglected. On the other hand, the spin- 
dependent lifetime, which represents the time for an in- 
version layer electron to irreversibly "fall" into the fer- 
romagnet, is in the picosecond range for oxide thickness' 
below 10 A, which approaches the intrinsic scattering 
time of the silicon inversion layer, tq « 1 ps^i The ratio 
of the lifetimes for the two spin channels t^/tj is roughly 
a factor of 2 for all oxide thicknesses; spin-dependent 
transport results, which can be harnassed for spintron- 
ics purposes. Using the following parameters for the 
Si/Si02/Fe interface, we calculate the scattering times 
to be tj = 3 ps and t\ = 6 ps. For silicon, the longi- 
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tudinal effective mass responsible for the confinement is 
m* n — 0.91m , the transverse effective mass responsi- 
ble for 2DEG transport is m* it — 0.19mo, the dielectric 
constant is e S i = 11.7, and the equilibrium density of the 
2DEG is assumed to be N a = 10 12 cm" 2 iS For silicon 
dioxide the effective mass is m£ = 0.3mo, the dielectric 
constant is eb = 3.9, and the barrier height (from the 
Fermi level) is Ub = 3.2 eVi^i For the ferromagnet, we 
use Uf m — 2.65 eV, exchange splitting A = 3.9 eV, effec- 
tive mass m* m 



TOor* and, for simplicity, a wavevector- 
independent imaginary part of the ferromagnetic energy 
7 | m = 1.1 eV and 7^ m = 0.8 eVM The electric field in 

the barrier is assumed to be E% = —10 MeV/cm. 

Having calculated the z-axis properties of the 
ferromagnetically-gated silicon field-effect transistor, we 
now examine the 2DEG transport in the indirection. Due 
to the electrical bias on the ferromagnet, the 2DEG cur- 
rent is not constant under the gate and will leak at a rate 
of 



d_ 

dx 



-N a (x) 



(7) 



Note that because the flow of electrons is opposite to the 
current, electrons leaking out of the 2DEG into the fer- 
romagnetic gate implies that a positive current flows into 
the 2DEG; the current will increase under the gate. Com- 
bined with Eq. and the appropriate boundary condi- 
tions (the 2DEG has its equilibrium density at the source 
and drain, N a (x s ) — N a (xd) = No/2), this results in 
a differential equation for the spatially-dependent, spin- 
dependent 2DEG density. The new scattering channel 
due to the interaction with the ferromagnet must also be 
included in the Drude conductivity, 



jV CT (x)e 2 ( 1 IT 1 



m 



si.t 



(8) 



The diffusion constant in two dimensions calculated using 
the above conductivity is 



D a (x) 



-F,a 



(x) 



<x)/k B T 



(9) 



where the spin-dependent Fermi level is ef a (x) = 
2irfi 2 N a {x) /m* it . In the following we assume the gate 
length in the indirection is approximately 1000 A, so that 
the source is at x s = and the drain is at x^ = 1000 A, 
and we assume that the gate extends 1 /jm in the in- 
direction. For concreteness, we will assume throughout 
this paper that the in-plane electric field is negative, 
E x < 0, such that, in the absence of leakage, a nega- 
tive current would flow through the 2DEG ( electrons 
flow from the source to the drain). In addition, we as- 
sume the confinement is homogeneous along the channel, 
so that Tf and are constant. 

In a silicon field-effect transistor with a thick gate ox- 
ide there is no leakage of 2DEG electrons because the 
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FIG. 2: (Top) A plot of the spin-dependent drain currents 
Jl(xd) (full line) and J^(xd) (dashed line) as a function of the 
in-plane field E x . (Bottom) A plot of the total drain current 
Jl(xd) + Ji{xd) (full line) and the spin current J] (xd) — Ji{xd) 
(dashed line) as a function of E x . At the in-plane fields 
marked A and C the drain current is 100% spin-polarized; 
at B, a pure spin current flows through the drain. 



wavefunction cannot penetrate into the gate. Neglect- 
ing the spin dependence for the time being, we examine 
the behavior of a silicon field-effect transistor with an 
ultra-thin gate oxide. Finite penetration of the 2DEG 
wavefunction through the oxide and into the gate causes 
a tunnel current to flow. The Fermi level in the gate is 
lower than the Fermi level in the inversion layer due to 
the gate bias; electrons that tunnel into the gate will in- 
elastically fall to the gate Fermi level and have no way to 
return to the 2DEG, decreasing the density in the 2DEG. 
The density will be near its equilibrium value No at the 
source and drain contacts, where electrons are injected 
into the 2DEG. The density will be lowest near the cen- 
ter of the gate, furthest from the source and drain. With 
no source-drain bias the density profile will be symmetric 
about the center of the gate; electrons will flow into the 
2DEG from both the source and the drain, such that the 
drain current is positive^ Application of a source drain 
bias will break this symmetry, causing a drift current to 
flow in addition to the diffusive current due the electron 
leakage. For negative in-plane field E x < 0, at some crit- 
ical field the drift current will exactly cancel the diffusive 
backflow at the drain; the net drain current will be zero. 
As the in-plane field is made more negative the drain cur- 
rent becomes negative, such that in the high source-drain 
bias regime the diffusive current is negligible. 

With a ferromagnetic gate, the above analysis still 
holds for each spin channel separately (provided that 
the spin-flip time is much longer than the transit time 
through the device). The lifetime for 2DEG electrons 
to fall into the gate is spin-dependent due to the spin- 
dependent coupling in the Hamiltonian, Eq. J3J). Im- 
portantly, the diffusive currents for spin f electrons and 
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FIG. 3: Plots of the spin-dependent currents Jf (full line) 
and J i (dashed line) under the ferromagnetic gate near field 
A (top), field B (middle), and field C (bottom), as labeled in 
Fig.H 



spin I electrons must also be different at zero source- 
drain bias, and hence the in-plane field required to ex- 
actly cancel the diffusive backflow at the drain contact 
will be different in the two spin channels. 

The spin-dependent drain currents Jj(xd) and J|(xd) 
are plotted as functions of the in-plane field \E X \ in the 
top panel of Fig. [21 The total current Jj(xd) + J|(xd) 
and the spin current Jj(xd) — Jj(xd) are plotted as func- 
tions of \E X \ in the bottom panel of Fig. El At very low 
source-drain bias the drain currents are positive in both 
spin channels, but have different magnitudes due to the 
different leakage rates for the two spin channels; because 
Tj < Tj, more spin f electrons leak than J, electrons. The 
larger gradient in the density in the spin f channel results 
in Jf(xd) > Ji(xd). At the in-plane field marked A, the 
drain current in the | spin channel is Jj(xd) — 0; the 
diffusive current is exactly cancelled by the drift current. 
The larger diffusive current in the spin f channel is still of 
higher magnitude than the drift current. Spin f electrons 
will flow into the 2DEG at the drain contact, resulting 
in a 100% spin-polarized current. At the in-plane field 
marked B in Fig. [21 the drain currents in the two spin 
channels are equal and opposite, J-r(xd) — — Jj.(xd)- No 
net charge current flows through the drain, but a pure 
spin current flows such that spin | electrons are trans- 
ferred into the 2DEG and spin J. electrons are transferred 
out of the 2DEG at the drain, i.e., the silicon inversion 
layer acts as a spin pump at the field marked B. At the 
field marked C the drain current in the spin j channel 
is negative, while in the spin | channel the drift and dif- 
fusive currents exactly cancel, Jj(xd) = 0. Spin J, elec- 
trons will flow out of the the 2DEG, resulting in a 100% 
spin-polarized drain current. As the magnitude of the 
in-plane field is increased further the electrical currents 
in both spin channels are negative; the spin J, channel 
carries more current because its conductivity is higher 



than that of the spin | channel. Electrons flow out of 
the 2DEG drain with some polarization. Note that the 
roles of the spin channels can be exchanged by revers- 
ing the gate magnetization, such that the drain contact 
of this ferromagnetically-gated silicon inversion layer can 
be used as a source for both fully spin-polarized charge 
currents and pure spin currents with either spin f or j. 
For positive in-plane electric fields E x > 0, the same phe- 
nomena happen at the source contact. 

The behavior of the spin-dependent currents through- 
out the 2DEG are plotted in Fig. [3| near the in-plane 
fields A (top panel), B (middle panel), and C (bottom 
panel). The leakage of electrons under the gate is evi- 
dent, as the current is not constant and increases under 
the gate. The plots show that in this region of source- 
drain bias the drift current is just beginning to overtake 
the diffusion current at the drain contact (the drift and 
diffusion currents flow in the same direction at the source 
contact). 

The spin effects present in the silicon inversion layer 
with a single ferromagnetic gate could be very attrac- 
tive for spintronics applications. We will now describe 
a simple all-electrical measurement to take advantage 
of the spin-dependent 2DEG transport. Recently we 
proposed a planar semiconductor spin-valve which relies 
upon the coupling of the inversion layer electrons to se- 
quential ferrromagnetsi22i2i The space between the fer- 
romagnets is filled with a paramagnetic metal to ensure 
homogeneous 2DEG confinement throughout the device. 
A diagram of the device is shown in Fig. \Rb) . The two 
ferromagnets are patterned with different aspect ratios, 
such that the magnetization of one ferromagnet can be 
reversed by an external field without affect the magneti- 
zation of the other ferromagnet. 

Assuming the two spin channels are decoupled 
throughout the device, let us examine the total 2DEG 
current at the drain contact for parallel ferromagnet mag- 
netizations, Jp(xd). One spin channel will see the spin 
t scattering time tj under both ferromagnets, and the 
other spin channel will see the spin J, scattering time t± 
under both ferromagnets; the total current J p (xd) will be 
the sum of the currents in the two spin channels. Switch- 
ing the magnetization of the second ferromagnet, one spin 
channel will see the spin f scattering time r-f under the 
first ferromagnet and the spin J, scattering time under 
the second ferromagnet, while the other spin channel will 
see the reverse. Because the density and current must be 
continuous in both spin channels, the total drain current 
for antiparallel magnetizations J ap (xd) will in general be 
different from J p (xd). The net result is a magnetore- 
sistance effect, in which the 2DEG current is dependent 
upon the relative orientation of the ferromagnet magne- 
tizations. We define the magnetoresistance ratio as 

MR = 100-| Jp(Xd) ~ J t p(Xd) | . (10) 

Jp(Xd) 

This quantity is plotted as the full line in Fig. 01 as a 
function of the in-plane electric field \E X \, for an infi- 
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where a is spin opposite to a. The magnetoresistance 
ratio is plotted for two values of the spin-flip time near 
the intrisic scattering time, r s f = 1.0 ps (clashed line) and 
r s f = 0.1 ps (dotted line) in Fig.0] Short spin-flip times 
drastically reduce the MR; the currents in the two spin 
channels are so strongly coupled that the total current 
through the 2DEG is only slightly altered when switching 
from parallel to antiparallel magnetizations. Although 
short spin relaxation times should not be a problem in 
siliconf2&2i it could be an issue in other semiconductors. 



10" 



|E J (V/cm) 



10 J 



FIG. 4: A plot of the magnetoresistance ratio (%) as a func- 
tion of the in-plane field \E X \ for the spin- valve with two fer- 
romagnetic gates, for different values of the spin-flip time r s f. 
For T s f = oo, the peak of the MR is limited only by the num- 
ber of points in the calculation. 



nite spin-flip time (in which case the peak of the MR is 
limited only by the number of data points in the plot). 
Both ferromagnets are assumed to be 1000 A long in the 
x-direction, with a 200 A gap between them. The mag- 
netoresistance ratio grows to values greater than 100% 
when Jp(xd) — ► 0, at which point a pure spin current 
would flow in the parallel configuration. 

The spin-flip time r B f couples the two spin channels 
together, driving them to the same density. Short spin 
flip times (r s f comparable to the scattering times To and 
t<j) will wash out any spin-dependent transport effects. 
For coupled spin channels, Eq. Q must be replaced with 



dx 



J a = 



N a (x) N a (x)-N s (x) 



(11) 



In conclusion, we have shown how the unavoidable is- 
sue of 2DEG carrier leakage in silicon field-effect tran- 
sistors with ultra-thin gate oxides can be put to good 
use for spintronics applications. By replacing the metal 
gate with a ferromagnet, the coupling between the sili- 
con inversion layer and the gate becomes spin-dependent, 
resulting in spin-dependent 2DEG transport. In the low 
source-drain bias regime, the drift current (first term on 
right side of Eq. QJ) and the diffusion current (second 
term on right side of Eq. 0J) are of the same order; the 
spin-dependent interplay between the two results in a va- 
riety of spin currents that flow through the contacts to 
such a device. We have shown how the ferromagnetically- 
gated silicon inversion layer can produce both fully spin- 
polarized charge currents and pure spin currents, without 
any net charge transfer. In addition, we have shown that 
if the gate is made up of two sequential ferromagnets, 
a magnetoresistance effect occurs in the 2DEG current 
dependent upon the relative orientation of the two fer- 
romagnet magnetizations. The magnetoresistance ratio 
can exceed 100% in this low source-drain bias regime. 
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